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Abstract. The competition between spin glass (SG), antiferromagnetism (AF) and Kondo effect is studied 
here in a model which consists of two Kondo sublattices with a gaussian random interaction between 
spins in differents sublattices with an antiferromagnetic mean Jo and standard deviation J. In the present 
approach there is no hopping of the conduction electrons between the sublattices and only spins in different 
sublattices can interact. The problem is formulated in the path integral formalism where the spin operators 
are expressed as bilinear combinations of Grassmann fields which can be solved at mean field level within 
the static approximation and the replica symmetry ansatz. The obtained phase diagram shows the sequence 
of phases SG, AF and Kondo state for increasing Kondo coupling. This sequence agrees qualitatively with 
experimental data of the Ce^Aui-xCoxSi^ compound. 
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1 Introduction 

It is recognized that there is a strong competition between 
the Kondo effect and the RKKY interaction in Kondo lat- 
tice systems . A transition from a magnetically ordered 
phase to a heavy fermion one, described by a Fermi-Liquid 
behaviour, has been observed in many Ce or Yb com- 
pounds and extensively studied from a theoretical point of 
view. There is a quantum critical point (QCP) at the tran- 
sition and non-Fermi liquid behaviors are also observed 
near the QCP 2 . The role of disorder has been studied by 
different approaches including a "Kondo disorder" model 
■ describing a broad distribution of Kondo temperatures (2j 
I or the extensive study of the so-called "quantum Grif- 
fiths" behaviour On the opposite, the transition from 
a metallic spin glass phase to a paramagnetic or Kondo 
phase has been studied recently by using the quantum ro- 
tor spin glass model and the existence of an anomalous 
behavior near T = has been observed in the transition 
between a metallic paramagnetic and a metallic spin glass. 

Particularly, spin glass (SG) and Kondo state have 
been found together in some Cerium alloys like CeNii-x 
Cux IEIj Ce2Aui-xCoxSi3 and in some disordered Ura- 
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nium alloys such as UCu^^xPdx H] or Ui^xLaxPd2Al^ 
In the first case, there is an antiferromagnetic (AF) 
phase for low contents of Ni. When the Ni doping is in- 
creased the phase diagram becomes more complex. For 
X < 0.8, the sequence of phases SG-ferromagnetism arises 
when the temperature is lowered and a Kondo state exists 
for X < 0.2. The alloys Ce2Aui-xCoxSi3 exhibit a phase 
diagram with a sequence of SG, AF and non magnetic 
Kondo phases with increasing the cobalt concentration at 
low temperature; in the high doping situation, the Neel 
temperature seems to tend to zero. In the two previously 
mentioned cases of disordered alloys 8 ,9 , the AF, SG and 
NFL phases have been obtained at low temperatures for 
different concentrations; in particular, the sequence AF- 
NFL-SG occurs with increasing x in UCu^-xPdx alloys [H] 
and the opposite sequence AF-SG-NFL with increasing x 
in Ui-xLaxPd2Al3 alloys 

Thus, the results previously mentioned evidence a quite 
complicated interplay among the Kondo effect and the 
RKKY interaction when disorder and frustration are pre- 
sent. Recently, a theoretical effort has been done to under- 
stand how a spin glass phase emerges in a Kondo lattice 
model with an intrasite exchange interaction and an in- 
tersite long range random interaction among the localized 
spins QSj^. The mentioned model has been extended to 
produce also a ferromagnetic order ^J. A particularity 
of our model compared to those discussed in Ref. 
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is that we do not consider quantum fluctuations by 
coupling only the Sz components of the localized spins in 
the spin glass term of the Hamiltonian, thus being unable 
to discuss a QCP at T = 0. One important remaining 
issue would be to obtain a description of the interplay be- 
tween antiferromagnetism, Kondo phase and spin glass. 
The first step in that direction would be to produce a 
theory at mean field level able to mimic some important 
aspects, for instance, the sequence of magnetic phases of 
Ce2Aui-xCoxSi3 in terms of a minimum set of energy 
scales related to the fundamental mechanisms present in 
the problem. In this work this set is Jq (the random inter- 
site average), J (the random intersite variance) and Jk 
(the strength of the Kondo coupling). In the following, 
we take a constant Jk and a Gaussian distribution of the 
intersite interaction with a random intersite interaction 
with average Jq and variance J. Thus, we will not take 
values of Jk and Tk given by a random distribution as in 
the work of Miranda et all W ■ 

In view of that, the aim of this paper is to present a 
mean field theory to study the interplay among disorder 
and antiferromagnetic ordering in a Kondo lattice provid- 
ing the necessary refinements in the previous theory [lUL 
The model studied is a two Kondo sublattice with 
an intrasite exchange interaction and an intersite random 
gaussian interaction only between localized spins in two 
different sublattices J^l- We consider here that there is 
no hopping of the conduction electrons between the two 
different sublattices; this assumption, which is used here 
to really simplify the calculations, does not modify the re- 
sults since the different AF or SG ordering are in fact es- 
sentially due to the localized spins. The possible AF or SG 
ordering is, therefore, entirely related to the coupling of lo- 
calized spins between the two sublattices. This fermionic 
problem is formulated by writing the spins operators as 
bilinear combinations of Grassmann fields and, through 
the static approximation and the replica formalism, the 
partition function is obtained, as already introduced to 
treat at mean field level the spin glass, Kondo effect and 
ferromagnetism in Refs. |TO] and ^1 . 

The Kondo effect and the RKKY interaction origi- 
nate from the same intrasite exchange interaction, but the 
necessity of considering an additional intersite exchange 
term has been already recognized and the full Hamil- 
tonian with both terms has been extensively used to de- 
scribe the Kondo lattice |14II15| . Moreover, the presently 
studied model has explicitly a random coupling term among 
localized spins in order to describe the spin glass situa- 
tion . The exchange terms cannot be considered as 
completely independent from each other and a relation- 
ship giving the intersite exchange integral varying as the 
square of the Kondo local exchange integral Jk has been 
introduced in order to mimic the {JkY dependence of 
the RKKY interaction ^Hl- In fact the relationship corre- 
sponds to an approximate representation of the intersite 
interaction, but we will use it in the last section of our pa- 
per in order to have a better agreement with experiment 
for some disordered Cerium alloys. Recently, the Doniach 
diagram has been revisited and the suppression of the AF 



and SG phases investigated, although these orders have 
been considered as two independent problems TF . 

This paper is structured as follows. In Section II the 
model is introduced and developed in order to get the 
free energy and the corresponding saddle point equations 
for the order parameters. In Section III, as mentioned in 
the previous paragraph, a relationship among Jq, J and 
Jk is introduced allowing to solve the order parameter 
equations and to build up a temperature versus Jk phase 
diagram showing the sequence of phases SG-AF-Kondo 
state. The conclusions are presented in the last section. 



2 General Formulation 

The model considered here consists in two Kondo sublat- 
tices A and B with a random coupling only between local- 
ized spins in distinct sublattices [12j. Furthermore, there 
is no hopping of conduction electrons between the sublat- 
tices as previously explained. The corresponding Hamil- 
tonian is given by: 

n-fiN^ H E XI ^ij4p,T^^j,P,<T + 5Z 4M,p 

p=A,B i,j CT=Ti i 



where i and j sums run over N sites of each sublattice. 
The intersite interaction J,., is assumed to be a random 
quantity following a gaussian distribution |12| 



^(-^-^^ jVei^'^n — 64j^ — ""j- 

An anti-ferromagnetic solution for the present choice of 
Jij can be found for Jq > 0. The case Jq < produces 
a complex phase diagram with spin glass, ferromagnetic, 
mixed phase (a spin glass with spontaneous magnetiza- 
tion) and a Kondo state as defined in the Ref. |1(J|. 

The spin variables present in the Eq. ^ are defined 
closely to Ref. (TUI modified for the case of two interacting 
distinct sub-lattices given as: sf^ — |di,p,xt ^ i^7p)^ ^ 

l' r ' ' 



fi,p,'\fi,pA fi, 



i,P,l 



^tp — fi,p,\fi,p,i — {^i,p)'^ andSfp 2 

/i,p,X where d[p,cr: dhp,<y {flp,a^ f^.p.'^) are the creation 
and destruction operators for conduction (localized) fermions. 

The partition function can be given by using Grassman 
variables i^i^p^aij) for the localized fermion and ^Pi^p^aij) 
for the conducting ones: 

•' p=A,B •' p=A,B 

where the action is A = + Ak + Asg with 

p=A.B i,j oJ 
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- E EE(4.pM)"'feTMV'. ,P:T('^) 

p=A,B ij 

+ ^lp,li^)^j,P,li^)] (4) 



and 



(5) 
(6) 



The energies a — b ~ ^^'^ A ~ 4 b ~ 4 
referred to the chemical potentials of the conduction and 
localized bands, respectively. 

The actions AgQ and Aj^ are given by 



-/ 



-/ 



^5g = /3EE^»^-^m(^)^^s(-^) 



(7) 



and 



A 



K 



N 



a=-\ip=A,B 



N 



N 



T.T.^lp>')^ 

i —1 UJ 



(8) 



with Matsubara's frequencies lu = (2m + l)7r and ly = 2m7r 
(rn = 0,±l,±2,...). 

The static approximation will be used in Eq. Q and 
Eq. ||SJ) to solve this problem at mean field level. In that 
spirit, the Kondo state at one particular sub-lattice s is 
caracterized by the complex order parameter Xp^a = jf 
J2i,uj (v*,p,ct{w)V'j,p,<t(w)) {p = A,B) which is introduced 
in the present theory through the identity 



N 



i.p.cri'^ 



27r 



exp < I 



E^p^' 



N 



^^;..-EE<p,-H^ 



UJ i—l 



(9) 



Its conjugate, A*^, can be also introduced by a simi- 
lar identity. From now on, it is assumed that \p a ~ Xp 

(a;,. « a;) [mil. 

Therefore, the partition function given in Eqs. ©-((SI), 
after using the integral representation for the delta func- 
tions, is 



Z = exp{-2NpJk{\\A? + |AbP)} 



(10) 



where 



^{stat)_ 



n -Dirp^p) I n T^i^ivp) 

=A,B P=A,B 
P [^0 



{stat) 
SG 



and 



Ak = Y,PJk 



^^^A.~ai^*^,A^)^P^AA^) + 
L lA 

^^\B,~ai^*jsAAVlBAA + 



>^*B,av)BAAi^]BAA 



(12) 



At this stage, the fluctuations in time and space are 
explicitly neglected. That means that in the sum over Mat- 
subara's frequencies, particularly for the spin part of ac- 
tion A^gQ*^ , only the term = is kept in Eq. Q . 

The conduction electrons can be integrated in Eq. (fTT)l 
to give 



^(stat) 



n V{rp^p)e^^^+^^" 

p=A,B 



where 



with 



^e//=EE^l(^ 



IjJ UJ .(T 



9. .{A 



(13) 



(14) 



{A 



[g%AAY'-FA{A 



[9%^b{AY''Fb{A 



f}^{A = { i^UAA ^IbAA ] 

^^(^^) ^ \ ^^^^AA ' 
[ ipj^BAA 

i^p(c.) = /3V||App^e^'^(^-^^)7fe(. 

fe 

where [5ijp('^)] ^ is given in Eq. JB)) and 



(15) 



ik{A 



1 



(icj - /3£o) - fifth 
The free energy is given by the replica method 

dF = fiJK{\\A\^ ^\\bA- 



(16) 
(17) 

(18) 



1 



lim 



- 1 



(19) 



where the configurational average ((...)) ^ is performed 
with the gaussian distribution given by Eq. ^ applied in 
Eq. which gives 

(11) {{A^"''AJ^,))) = I D{rAA^"{{z,,)) (20) 
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where 



{(Zso)) 



G//CA 




Therefore, the resulting averaged partition function 
can be hnearized by using the usual Hubbard-Stratonovich 
transformation and introducting auxiliary fields of the spin 
glass part 12 . The details of this procedure are given in 
the Appendix. 

The free energy can be found using the averaged parti- 
tion function (see Appendix) in Eq. (fT^ . Thus, one gets, 

J2 

+ P^-^i^AQB) ~ PJomArriB 

I l^oo ^ poo 

- }^^o 2N-n y_ n n y_ 



X exp 



where the 



^Indet 



(22) 



G.\u}\hf ) in Eq. ^ is given by 

1 J 



{icu- 







{iu;~p4~ah-s)S,^,-FBiu;) 



(23) 



with Fp defined in Eq. H17|l and hfj^ and hfg being ran- 
dom gaussian fields (see Eq. H34|l from Apendix) applied 
to the sites of one sublattice {A or B) which depends on 
the parameters of the other sub-lattice. At this point, the 
static susceptibility Xp can be introducted related with 
replica symmetry diagonal order parameter (jp = Xp + Qp 
where Xp = Xp/P- 

The central issue here is to adopt the proper decou- 
pling approximation which allows one to calculate the ma- 
r 1 ~^ 

trix G . . (w) nUTT] . The elements given in Eq. are 

referred to the original sublattice A (B) where to each site 
i there is a random gaussian field /i"^ (^"b) a-pplied. The 
decoupling procedure is to consider these random fields 
^7 A (.^?b) ^ apphed in two fictitious Kondo sub-lattices. 
Therefore, in each site fi of the new sublattice, the applied 
field (^"b) constant. That is equivalent to replace 



by [£,MKp) 



the Green's function G..{oj\hf ) 
L — 

in Eq. (|23|) . It is, therefore, possible to go to the reciprocal 
space, where it is assumed a constant band p{e) = for 
—D < £ < D ior the conduction electrons in each sub- 
lattice. The sum over the Matzubara frequencies in Eq. 
H22(l can be done following a standard procedure. Thus, 



the free energy is found to be 



PF = PJK{\\Ar + |Asr) + ^-l^XAXB + ^-^iXAQB + 

XBqA) —mAmB - x / D^^^ / D^j^ x 



In 



n 

p=A.B-'^°^ 



(24) 



with 



1 r'" 



, , x + hp 
cosh ' 



cosh 



' X — hr 



(25) 



The saddle point equations for the order parameters 
%: iTT-pj Xp and |Ap| 10 follow from Eqs. (|24|) and H25|l . 



3 Results 

The numerical solutions of equations for the order param- 
eters Qp, Xp, nip {p — A, B) and |Ap| allow one to build 
up a phase diagram of temperature T versus the set of 
relevant energy scales Jq, J and defined in Section II. 
The thermodynamic phases are identified as: (a) an anti- 
ferromagnetic (AF) phase corresponding to ruA = —tub 
or 7713 = (see Eg. I^ - IH^ ): (b) a spin glass phase for 

7^ and g_B 7^ (ga ^ 0); (c) a Kondo state given by 
|Aa| ^Oand |Ab| ^ 0. 

The phase diagrams obtained by the present calcula- 
tions are derived for several sets of parameters Jk/J and 
Jo /J considered firstly as independent from each other. 
But, as it was discussed in the introduction, the intrasite 
and intersite exchange terms cannot be considered as com- 
pletly independent from each other and a J|-dependence 
of Jo was introduced in order to mimic the RKKY inter- 
action JHl- Thus, a phase diagram is shown in Fig. |^ by 
keeping J constant, taking then: 



Jo 
J 



(26) 



and by choosing finally a — 0.0051. The factor a was 
chosen here to have the AF phase starting at Jk/J ~ 12. 

The numerical solutions of the order parameters yield 
the sequence of phases SG-AF-Kondo shown in Fig. (QJ. 
The Kondo temperature Tk is decreasing for decreasing 
values of Jk/J down to J^/ J Ki 15. ^From this point on- 
wards there is a phase transition leading to antiferromag- 
netic order with the Neel temperature T/v monotonically 
decreasing with Jk/J- For still smaller values oi Jk/J (and 
therefore smaller values of Jq / J) the spin glass behaviour 
becomes dominant. 

Thus, in this mean field theory it is possible to identify 
three sorts of distinct regimes depending on the strength 
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- Jo/J=0.0051(Jk/J)^ 
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Kondo 






A AF \ 




sg 

1 1---'"' 1 
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10 



20 



(Jk/J) 



Fig. 1. The phase diagram T/J versus Jk/J showing the se- 
quences of phases SG (Spin Glass), AF (Antiferromagnetism) 
and a Kondo state. The variance J is kept constant. The dotted 
hne means the "pure" Kondo temperature. 



of the Kondo coupling J^- The first one (for high Jj.) is 
the complete screening of the localized magnetic moments 
due to the Kondo effect, which is caracterized by the or- 
der parameter |Ap| related to the formation of d-f sin- 
glet throughout the whole two sublattices. In the second 
one (for decreasing J^) , the magnetic moments of one 
particular sublattice, for instance sublattice A, survive to 
the screening process (the complete unscreening means 
|A^| ~ and |As| ~ 0) and start to be in an antiparallel 
alignment with an internal field /i^, which depends on the 
magnetization and susceptibility of the sublattice B, until 
this internal field is spread out to the entire sublattice pro- 
ducing the AF order. In the last regime ( for J^/J < 12 or 
Jo < 0.734 J), the effects of randomness begin to be dom- 
inant. In the effective field hp [p — A^B), the replica spin 

glass order parameter q^i {p = B, A) component starts to 
be non-null indicating the non-trivial ergodicity breaking 
leading to a spin glass phase at the transition temperature 
Tf. Finally, let us remark that we never obtain theoreti- 
cally mixed phases where there is the coexistence of the 
Kondo phase with the other phases. This is due to the 
mean field approximation used here, as already observed 

The experimental phase diagram of the alloy Ce2Aui-x 
C Ox Sis |7j can be addressed if the Jk coupHng is associ- 
ated with the content of Co. The obtained result shown in 
Fig. ^ displays the same sequence of phases at low tem- 
perature as the experimental one. In brief, the increasing 
Jk favours the transition from a SG phase to a AF order- 
ing and then to the screening of the localized moments. 
Thus, this mean field description is able to account for 
experimental aspects of Ce2Aui^xCoxSi3 at low temper- 
atures. 

There is some disagreement related to the location 
of the AF line transition. The experimental behavior |7j 
shows the Neel temperature Tjv decreasing apparently to- 
wards a quantum critical point (QCP) with increasing Co 
doping. Howewer, recently the spin flipping has been sim- 



ulated with the presence of a transversal field in the x- 
direction to study a QCP in fermionic spin glass ^Jj. The 
same approach has been extended to investigate the spin 
glass freezing in a Kondo lattice Hence, this method 
could also be used in the present model in order to clarify 
the role of the QCP in the interplay between spin glass 
and antiferromagnetism in a Kondo lattice. That will be 
subject for future investigations. 

In conclusion, we have studied in detail the competi- 
tion between SG, AF and Kondo phases in a mean field ap- 
proach of Kondo-lattice disordered alloys. Using the pecu- 
liar relationship (|26|l between the different parameters of 
the model, we have obtained a phase diagram showing the 
sequence SG-AF-Kondo phases in good agreement with 
the experimental phase diagram of the disordered alloys 
Ce2Aui-xCoxSi3, if we assume that Jk increases with 
increasing Cobalt concentration. Other phase sequences, 
involving also SG, AF and Kondo phases, obtained in dis- 
ordered Uranium alloys, are not accounted for, at present, 
by our model. But the choice of parameters and their re- 
lationship with the varying concentration in these alloys 
are really delicate questions which are not answered at 
present. Moreover, in such Uranium systems, the Kondo 
phases have clearly a Non-Fermi liquid behavior which is 
not described here. Finally, the present model can account 
for the presence of SG, AF and Kondo phases, but further 
work is necessary in order to obtain a better agreement 
with experimental data for Cerium and Uranium disor- 
dered alloys. 
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Appendix 



The averaged partition function can be given by |12| : 



p=A,B a 

X ' n n 

p=A,B a/3 



A D r^, * 



2tt ^^-p 



a/3 



27r 



rfQf /n 



27r 



dM^ 



X exp \ —N 



--In AiM^,Qf) 



Q/3 



p=A,B 
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/3Jo 



D=A.B 



(27) 



where the A{M!^, Q'^^) is defined as 



A = J D{ri^)e^plY.^H 



+ iPJo E 

a. 

+ /3^oE 



a/3 



1 



5'..('^) 



M3" 



a/3 



E/ "^i^^^a + E/ ^j'B.a 

j 

'y^,{^i.A.a^i.A.0)Q A ~^ 

j 

E '^J\^,a'5'i',A,/3 + 
E'^I^S,a'S'|,B,/3 




(28) 



/From Eqs. (|77|l and if^ . one finds the saddle point 
solution for these auxiliary fields I12j : 



M!; 



1^(T.^Ia,.^11^1b.o) = '^^rn'i (29) 



N 



K = ^ ( E ^^P." / = ; P = ^, s (30) 



~ \ E ^i,A,a^i,A,P + E ^j,B,a^: 



4i 



'j,S,/3 



n ^^^^ n / exp [/i^s5r,B,a] 

/OO p OC 

U^^^^Ii Dzr^^exp[h1^^SUa] (33) 



(31) 8 



~ {^^^^i,P,a^i,p,Pl — '''Ip'^' P — ^7 B (32) 

It has been assumed, within the replica symmetry ansatz 
for the auxiliary fields (see Eqs. 129() - (|32() '). that q^^ = q^, 
Qp^ — Qp and (7p" = (jp (analogously = and 
rrip = nip). The sum over replica index produces again 
quadratic forms which can be linearized by new auxiliary 
fields. Thus, for yl(m3,g3) one has 



j D{r^)e^A E ^U(^) 



where the random field hf^, introduced in the previous 
equation, is defined as 

htp = PJ {^^^,p + ^J^Qp' -qp')ztp) 

-(3Jomp' (p^p) (34) 

with Sl^^^ = IJ2 '^{^tpJ*i^mpA^) and Dx = 

(e-"'/V\/2^) dx. 

The functional integral in Eq. (|33|l can be found fol- 
lowing standard procedure |19 |. 
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